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Higher-order diagrams required for radiative corrections to mixed electroweak and QCD processes at the LHC 
and anticipated future colliders will require numerically stable representations of the associated Feynman dia- 
grams. The hypergeometric representation supplies an analytic framework that is useful for deriving such stable 
representations. We discuss the reduction of Feynman diagrams to master integrals, and compare intcgration- 
by-parts methods to differential reduction of hypergeometric functions. We describe the problem of constructing 
higher-order terms in the epsilon expansion, and characterize the functions generated in such expansions. 
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1. Introduction 

A variety of approaches are known for evaluating Feynman diagrams analytically [l[ and constructing the 
e expansion of dimensionally-regularized integrals. One-loop diagrams are expressible in terms of 

generalized hypergeometric functions of the form p +\F p , Appell functions Fi, F 2 , F 3 , F 4 , or Horn-type 0] 
multivariable functions. [8| Analytical techniques for evaluating the finite part of one-loop diagrams have long 
been known, Q but higher order diagrams require more terms in the s expansion. These facts motivate us to 
seek a way to obtain the all-order e expansion of Horn-type hypergeometric functions in general. 

The first such construction of the all-order e expansion was obtained for the Gauss hypergeometric function 
with special, physically-motivated, values of the parameters. 10] The first systematic algorithm for construct- 
ing the e expansion of multi- variable hypergeometric functions was described in Ref. [12j |. using Goncharov's 
poly logarithms [ill ] in the case with integer values of parameters. However, diagrams with massive propaga- 
tors lead to hypergeometric functions with rational values of parameters. 11311 This case generates multiple 
(inverse) binomial sums which were first systematically investigated in Refs. [14 - lla ]. Ref. [T7] generalized the 
technique of Ref. 12] to the case of rational parameters, but the method was limited to the "zero-balanced" 
case, and applied to p +\F p and Fi, but not F 2 or F 3 . Some statements about the structure of the coefficients 
of the e-expansion for p +\F p with one non-balanced rational parameter contradicted explicit calculations of the 
first few coefficients of the e expansion. [l(| [l8| In particular, the multiple sum Y^jLi -^b'C/)] -1 Y^k=\ l ^ 2 
was evaluated analytically in Ref. [l6[ (see Eq. (C.22) there), and it was found that the e expansion of the 
hypergeometric function 



P+iFp 



h + h + fie, h 



b\E , Kj + djE 



Ci£ 



(1) 



with integer values of Ij,Kj,Lj leads, at weight 3, to classical polylogarithms of argument jqr^f- (This was 
later confirmed in Ref. [l9j.) Furthermore, Ref. [HI analyzed the more general case 



r , h + | + bie , I 2 + \ + b 2 e , Kj + aj£ 

p+1 P \h + k + fie, h + h + f 2 e, h + h + he, U + c, : 



(2) 



via the sum Y^JLi z ^\j{ 2 j)] 1 Y^l—i 1 k 3 , and was again found to differ from the prediction of Ref. fl7l ]. 

The inadequacies of existing algorithms and a desire to explore the results of Ref. 0] for one-loop diagrams 
motivated the development of new technologies (20l - |25j and the extension of previous ones [l6[ for all-order 
expansions, [26[ with the goals of 

(i) independently verifying the results of Ref. [TtJ , 

(ii) constructing the analytical coefficients of the e expansion of hypergeometric functions of several variables. 
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Recently, new _generalizations of the algorithm of Ref. 12j have been found using cyclotomic harmonic sums [27] 
and Z-sums, [28] allowing the e expansion to be constructed in the cases of Eqs. ((T|) and ((2J. Their applicability 
to multivariable hypergeometric functions is not yet clear. 

The higher terms in the e expansion of one-loop diagrams can be constructed Q without applying algorithms 
for expanding hypergeometric functions. [29|-[32j However, results derived via the hypergeometric representation 
are more compact and often can be expressed in terms of simpler functions, e.g. Nielsen's polylogarithms instead 
of Goncharov's. This can be seen, for example, by comparing the results presented in Refs. [29|, [3(| to those in 
Refs. [lol [lH [33j. Some remarkable results on hexagon diagrams were presented in Ref. (32| . but this was done 
via an experimental mathematical analysis [34[ specialized to this case. Methods following from the properties 
of hypergeometric functions hold the promise of being applicable to wide classes of Feynman diagrams. In 
addition, the e expansion of hypergeometric functions has interesting applications beyond Feynman diagrams, 
in a broader mathematical context. [35[ 

In this paper, we will concentrate on Horn-type hypergeometric functions, a very general class which encom- 
passes all of those mentioned earlier, and which we conjecture [IB] to be general enough to express all Feynman 
diagrams. Section [2] describes the differential reduction of this class of hypergeometric functions. Section [3] 
describes the differential reduction in application to Feynman diagrams and explain the difference in count- 
ing master integrals via our differential technologies [37| versus integration- by-parts techniques (38| . Section [4] 
describes our approach to constructing the e expansion, and provides some examples. 

2. Differential Reduction of Horn-type Hypergeometric Functions 

Let us begin by considering the generalized hypergeometric functions p +iF p (a; 6; z) defined by a series about 
z = as 

^ p ( , M _^____, (3) 
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where {a)u = T(a + k)/T(a) is the Pochhammer symbol, and a — (<zi, • • • , a p ), b — (b\, ■ ■ ■ , b q ) are called the 
upper and lower parameters, respectively. This function satisfies a differential equation of order p + 1, so that 
there are p + 1 independent solutions for unexceptional values of the parameters. 

"Differential reduction" is an algorithm by which any function p+ \F p {a + m; b + k; z) with integer parameter 
lists m, k may be expressed in terms of a set of p functions in which the values of these arguments are shifted 
by lists of integers (l)j and (r)j, respectively, with j = 1, . . .p. This reduction has the form 

p 

Q p+ i(a,b,z) p+1 F p (a + m;b + k;z) = ^Q j (a,b,z) p +iF p (a + lj,b + rj; z) (4) 

3=0 

where the Qj are a set of polynomials in ai,bi, and z. The resulting expression can be converted into derivatives 
of the unshifted hypergeometric function, so that the reduction takes the form 39] 

p 

S(a, b, z) p+ iF p (a + m; b + k; z) = Rj(a, b, z) 9 j p+ iF p (a; b; z) , (5) 

3=0 

where 9 — zd/dz and S and Ri are other polynomials in ai,bi, and z. 

For particular values of the parameters, an algebraic solution of the corresponding differential equations can 
be found, so that the differential reduction simplifies. For example, when some of the upper parameters are 
integers 7j, we have [4(| 

p-i _ 

P(a, b, z) p+1 F p (I, a + m; b + fc; z) = ^ Rj(a, b, z) 9 J p+1 F p (I, a; b; z) + R p (a, b, z) , (6) 

3=0 

with algebraic functions R p (a, b, z). In this case, the Eq. ([4]) has the form 

p-i _ 

Q p+ x(a, b, z) p+1 F p (I, a + m; b + k\ z) = 2J Qj(a, b, z) p+1 F p (J, a + lj;b + r d \z) + R(a; b; z) , (7) 

3=0 
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where I, J are lists of integer parameters. 

This reduction can be generalized to all multivariable Horn-type hypergeometric functions (39] having L 
independent solutions: 

(i) For unexceptional values of parameters (irreducible monodromy group), linear differential relations can 
be found among L + 1 functions with parameters differing by integers. 

(ii) When the monodromy group is reducible, the number of functions entering into these differential relations 
is correspondingly reduced. 

In the reducible case, a basis can currently be constructed only when the differential reduction has been carried 
out explicitly. 



3. Counting the Number of Basis Elements 
3.1. Preliminary Considerations 

Consider a function f(z) satisfying a homogeneous differential equation of order k with polynomial coefficients: 



$>(*)(-£) m=o. (8) 



1=0 



It is easily shown that the function H(z) = z a f(z) satisfies a differential equation of the same order as f(z). 
Consider a function H(z) which is a linear combination of functions fj(z) satisfying homogeneous differential 
equations of order with rational coefficients, 

m 

H(z) = J2^x ak fk(z), (9) 
fc=i 

(|)"aw+E^)(s) , ^ = - (10) 

where ctk are rational numbers and r^ are integers. Then it can be shown that the function H(z) satisfies a 
differential equation with rational coefficients of order not less than m x r^. 

We can take the fk to be Horn-type hypergeometric functions (multi-variable, in general) whose parameters 
are linear combinations of the parameters of some Horn-type function H{J; z), which plays the role of H(z) 
above. Then, in accordance with Sec. [2] there are linear differential relations between any m x rfe + 1 of the 
functions H(J + m; z). This means that these functions can be expressed in terms of a set of basis functions 
shifted by integer values of the parameters, 

The number L of elements on the r.h.s. of this relation is equal to the number of solutions of the corresponding 
differential equation. 



3.2. Application to Feynman Diagrams 

As an application of the reduction in Sec. [21 a series of publications [40M43l | has analyzed the Feynman 
diagrams &(n,j; z) having a Mellin- Barnes representation of the form 

$(n,j;z) = — dt(KzYT(~t)F(t) , (12) 

where j is the list of powers of the propagators in the Feynman diagram, n is the dimension of space-time, n is 
a a constant, and F(t) has the structure 

n r(A l + t)r(c k -t) 
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where Aj, Bj, Ck, Di are linear combinations of n and j of the form A4 = ao n +X)fc a kjk , etc., and the dimensions 
of the lists A, B, C, D of these linear combinations satisfy the condition 

dim A + dim D - dim B - dim C = 1 . (14) 

Closing the integration contour gives 

q 

<f>(n,j;z) = ^z l °C la {n,j) p+1 F p (A a ; B a ; kz) , (15) 
0=1 

where q is an integer, l a , A a , B a are linear combinations of n and j with rational coefficients, z / 1, and the 
coefficients Ci a are products of T- functions with arguments depending only on n and j. 

In accordance with the preliminary remarks in this section, the maximal number of basis elements for Eq. (|15|l 
is equal to q x p. This number should coincide with the number of master integrals found by other means, in 
particular, by integration- by-parts (IBP) relations [381 ] . However, some of the parameters of the hypergeometric 
functions on the r.h.s. of Eq. (|15p may be exceptional. Typically, this occurs when the upper parameters are 
integers, or the difference between an upper and lower parameter is a positive integer. [4fJ In these cases, the 
actual number of basis elements will be less than the typical number. (See also the discussion in Ref. (42j.) 

The number of nontrivial basis elements for a hypergeometric function of one variable may be defined to be 
the highest power of the differential operator 6 = zd/dz in its differential reduction: 



p+1 F p {A;B;z) = ^ 



J2 p i( z ) ° l s+iF s (A - h, B - I 2 , z) + R s {z) 



.1=0 



(16) 



where ii, J2 are lists of integers and Pi(z), R(z) are rational functions. On the other hand, applying the IBP 
relations to Eq. (TTBl would lead to an expression of the form 

h 

®(n,j;z) =J2 B k( n J: z )^k(n;z) , (17) 
fc=i 

where some of the master integrals $fc are expressible solely in terms of r-functions. [44| 

Analyzing the relation between the reductions ((16)) and (IT71) leads to the following criteria for determining 
the number of master integrals: 1 

(i) Each term in the hypergeometric representation of a Feynman diagram Eq. (j!5[) has the same number L 
of nontrivial basis elements (up to rational functions). 

(ii) The number of master integrals following from the IBP relations which are not expressible solely in terms 
of r-functions [44j is equal to the number L of nontrivial basis elements. 

3.3. Examples: Vertex-Type Diagrams 

In this section, we consider two examples of the differential reduction of one-loop vertex-type diagrams. 
Example I: Let us consider the one-loop vertex (C3 in the notation of Ref. |40j : see Eq. (52) there for details): 



d n k 



[(fc - pi) 2 - m 2 p'i [(Jfc + p 2 ) 2 - m 2 }iz (k 2 Y 



.-i-n_n/2, m h^-fa F (jig + ° - f ) r (j - a) ( j 12 + a -%,ji,j 2 ,z-a 



71 y~ m ) 2 r /~ ^.wm 4 3 n J' 12 212+1 

1 Ul2)l (j) \ 2' 2 ' 2~ 



(Pi -V2)' 

4m 2 



(18) 



where j'12 = ji + J2, v = Sfcii a k, and the "dressed" massless propagator is 



1 _[ff r(§-^) ] ^-"W^r^-fg) 1 



(fc 2 ) CT U=i r(CTfc) J r(f(9 + i)-^) (fc 2 r-f" 



1 See Ref. [40ll for details. We thank Bas Tausk for bringing to our attention the fact that the diagram F2 in that paper was 
analyzed in detail in Refs. [4S| (Eqs. (74), (75)) and [49| (Eqs. (124), (130)). The results of these analyses are in agreement. 
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where q is the number of massless loops, (q + 1) is the number of massless lines, and q — corresponds to a 

massless line without additional internal loops. (Cg ^ is a true one-loop vertex.) 

With the redefinition a — > a — '^q, the hypergeometric function in Eq. (fT5| takes the form 



4^3 



j + a-%{q + l),j u j2,%(q+l)-<T 

n £ j + 1 
2 ' 2 ' 2 



(Pi -P2? 

4m 2 



(20) 



In accordance with the differential reduction algorithm, this function may be expressed in terms of a 3P2 function 
with one unit upper parameter, and its first derivative, 



{l,6}x 3 F 2 



l,5(« + l) + Ii,/ 2 -f(«+l) 

n 1 
2 ' 2 



n 1 1 t 
9 > 2 + J 3 



(21) 



together with rational functions of z, for integer values I\, J 2 , ^3 ■ The short- hand notation (1,0) stands for a 
combination Pi(z) + P2(z)9, with rational functions P$. (See Eqs. (17) and (20) in Ref. [4y| .) According to 
the criteria in Sec. 13.21 there are two master integrals for this diagram which are not expressible in terms of T 
functions. For q = 1, the standard approach based on IBP relations, [H| yields two two- loop vertex master 
integrals. (See Eqs. (3) and (9) in Ref. [5(|.) These master integrals are relevant for the massless fermion 
contribution to Higgs production and decay. [5l[ It is interesting to note that at the one-loop level, 
has one master integral of the vertex type and one of the propagator type, which are again equivalent to Eq. (|21j) . 



Example II: Let us consider one-loop vertex diagram C\ 



, defined as 



d n k 



[(k - pi) 2 ] CTl [{k + p 2 ) 2 ] CT2 (k 2 - m 2 )P 



(22) 



=pI=o 



where «n, er 2 are defined as in Eq. (|T9|) . with Oj = Y^k^i a i k f° r 3 — -*-! 2- The case q% = 0, g 2 = 1 corresponds 
to Eq. (173) in Ref. The hypergeometric representation for this diagram is 53] 



(11,12,) 



(<Ti,<T2,p)=e-V(-m')i 



2\-^ — p—o 1 -o 2 



r(p+ C r 1 + ( 7 2 -5)r(|-a 1 -a 2 ) 



r(f)r(p) 



J-F2 



P+CTx + O-2- 2,CTi,CT 2 
-, l + cri + a 2 --5- 



(pi -P2) 2 \ 2 



X3P2 



U — (Ti — <t 2 , 1 — Or — cr 2 + 1 



r (f ^^1) r (f r (<x 1+0 - 2 -f ) 

r (7i-(7i-(T 2 )r((Ti)r((T 2 ) 
(pi -p 2 ) 2 \ I 



(Pi -P2) 1 



(23) 



In accordance with the differential reduction algorithm, the first 3P2 function in Eq. (|24p can be expressed in 
terms of 2 Pi, and the second one can be expressed in terms of 3P2 with one unit upper parameter, namely 



{1,6} x 2F1 



0\,02 

n 
2 ' 



{1,0} X 3^2 



in _ n _ 
J-) 2 — CTl, 2" 2 

n-CTl-(T2, f-CTl- CT2 + 1 



(24) 



together with rational functions of z. We note that each hypergeometric function in Eq.[53]has the same number 
of nontrivial basis elements, in agreement with criterion (i) in section 13.21 When one (or both) of the o~j are 
integers (qj = 0), further reduction is possible to 



1,0-2 

n 
2 • 



2 Pi 



1, 3-^2 

n — o\ — o 2 



(25) 



In accordance with our criteria, there are two nontrivial master integrals for non-integer oj, and one nontrivial 
master integral when one of the o~j is an integer. The last result in agreement with Ref. [52j . 
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3.4. Counting the Non-Trivial Master Integrals 

When analyzing various diagrams, one example was found [37l | where our criteria allow us to predict and prove 
an additional relation between master integrals. It is a two-loop sunset diagram relevant for the evaluation of 
0(aa s ) relation between the pole and MS mass of the top-quark in the Standard Model, [54j as well as in the 
gaugeless-limit model. [55j To illustrate this, let us consider the two-loop self-energy diagram 



V 1200 (p, a,a,P,m 2 ,M 2 ) 



d n {k 1 k 2 



[(fc 1 -p) 2 -m 2 ]"[(fc 1 -fc 2 ) 2 - M 2 }°[k 2 } a [k 2 



2 V j 



(26) 



which is generated in an intermediate step of the calculations in Ref. [55]. The a = case corresponds to a 
diagram considered in Ref. |37[. The Mellin-Barncs representation of the integral (|26p is 



V 120 o(p,<y,a,l3,m 2 ,M 2 ) = [i^V 



m-a) 



-(-m 2 ) 



2\n—a—fi—a—p 



ds 



M 



J2 



T(a)T(p)T(a) 

T(-s)T(^-(7-s)T{2n-2a-2l3-2a-p-2s)T(a+f3+a+p-n+s)T(a+a-^+s) 
T(n-a-a-s)T(^-a-f3-a-p-s) 



(27) 



After integration, we obtain 



Vi2oo(p, o-, a, P, m 2 ,M 2 ) = [i 



l-n 



7T2 



r(|-a) 



771 

~M 2 ~ 



r( a )r(a)r(f 

~ ff - p T(a)T (a+(T-|) r {P+p-f) T{n-2(3-p) 



f_ jyr2\n-a-f3-<T-p 



T(p)T(n-p-p) 
r {a+13+a + p-n) V (a + /3 + p-f ) T (j-p-p) 

ne+p) 



4-F3 



n+l—p _ 

f,n-/3-^| + l-/3-p 



■ n n — p n 



r , a + f3 + <T+p-n,a + f3 + p- 7 j , f , £±± 



Am 2 
M 2 " 



M 2 



(28) 



The differential reduction of the hypergeometric functions in Eq. (|28|) can be expressed as 



4^3 



4^3 



a, a+c- 
" ,n- 



_ n n — p 
' 2 ' 2 



2 ' ^ f ' 2 

a + /3 + cr+p — rt,a + /3 + p 
f,/3 + p,l+/3+p- 



P-P 



1- 



-> (1,0) x 3F2 



n £ p+1 
2 ' 2 ' 2 



1 j n n _ 

-L 1 * 1 9 5 9 



n n 
2 ' 2 



■3 + 

-/4 



(M) X 3*2 



/i-n, 



-/ 2 ,1 



(29) 



In accordance with our criteria, there should exist two algebraically independent master integrals and some 
integrals expressible in terms of rational functions. 

Diagram (|2"6"f is algebraically reducible via IBP relations [47[ to four new diagrams with (i) a = 0, (ii) p = 0, 
(iii) a = 0, (iv) P = 0. The <r = diagram vanishes in framework of dimensional regularization. The diagrams 
with P = and a = are irreducible, and treated as two independent master integrals. The diagram with p = 
is a two-loop sunset, and in accordance with the IBP algorithm of Ref. |47| . it has three master integrals of the 
same topology plus a product of two one-loop bubbles (a = P = 0). In this way, the classical IBP relations 
applied to diagram (|26[) give six master integrals, without any information about its algebraic structure. The 
hypergeometric representation for sunset diagrams (|2"5|) and the differential reduction algorithm (JSJ allow us to 
find algebraic relations between the three sunset-type master integrals. (See Eq.(9) in Ref. 37[.) The on-shell 
case of a new relation (Eq.(10) in Ref. [37]]) enters in the package ON-SHEL L2 |56l . and was postulated via a 
study of the higher-order coefficients of the e-expansion in Ref. [13] • In Ref. |58|. it was pointed out that the 
last on-shell relation can be treated purely diagrammatically. In this case, the new diagram does not follow 
from the original sunset by the contraction of any lines. Surprisingly, this new relation is not reproduced by 
the IBP algorithm of Ref. [47j] or by the automated programs AIR [59| and FIRE [goj|. This may be an artifact 
of the "topological" nature of the solution of IBP relations (systematic contractions of one line in a diagram). 
The effect that the solution of IBP relations does not recognize the algebraic relation between master integrals 
is also seen for phase-space integrals: some of the master integrals collected in Ref. [6l| are also algebraically 
related to each other. 
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4. Construction of the e Expansion 

The evaluation of multi-loop radiative corrections requires higher-order terms in the e expansion of lower- 
order diagrams. The i nteg ral representation (53|, the series representation [l2], HH-Qjl Hfij or the differential 
equations approach [2CJ, l23l425l | may be used to construct the e expansion of hypergeometric functions. We will 
focus on the differential equations approach. As is well-known, hypergeometric functions satisfy a differential 
equation 



p F p ^(A;B;z)=0. 



(30) 



Due to the analyticity of the hypergeometric function p F p -i(A; B; z) with respect to its parameters Ai,B k , 
the differential equation for the coefficient functions u k of the Laurent expansion can be derived directly from 
Eq. pop without any reference to the series or integral representation. This was the main idea of the approach 
developed in Refs. ^[H, Hi- 
lt is convenient to introduce a new parametrization, Ai — > Ai + ct^e, Bj — > Bi + bie , where e is small, so that 
the Laurent expansion takes the form of an "e expansion," 

oo 

u){z) = p F p _i(A + as; B + be; z) = p F p - X {A) B, z) + ^e k uj k {z) . 

k=\ 

We can rewrite Eq. (f3"U|) as a system of first-order differential equations (called the "Pfaff form") and expand 
all terms in powers of e: 



p-i 



d4> (l) (z,e) = J2 ( z > ( z > £ ) dz 

3=0 



(31) 



where <j)^(z) — hi(z)9<j)^ ^ (z), i = 1, • • • ,p— 1 and = uj(z), with arbitrary functions hi(z). The result can 
be written in the form 



(32) 



where Pk{z) and R^( m \z) are polynomials. For a restricted set of parameters, the expanded system of dif- 



ferential equations in Eq. (|32[) has triangular form and can be integrated iteratively in terms of Goncharov's 
poly logarithms . 2 

For illustration, let us consider the factorization of a differential operator expanded in powers of e: 



d (p) 



1 



Ai + die) — 9IL p k Z\ {0 + B k -l + b k e) 



j2eW^\A,B,a,b,z), 

3=0 

where 9 — zd/dz, the upper index gives the order of the differential operator, Dp " 1 = Tl k=1 a k , and 



(33) 



D, 



(v) 



nLi 



-Ai)--eil p k Z\ (9 + B k -l) 
p-1 



a-^+E^Ec^- 1 ) ^- X +E 



fe=i 



3=1 



3 = 1 



(34) 



i=0 



with polynomials pj p \r\, ■ ■ ■ ,r p ) defined via the relations 

TJ (z + r k ) = jr pW ( ri , • • • , r^E £ P p W (r)z^ £ P^\r) 



yP-3 



k=l 



3=0 



3=0 



3=0 



(35) 



2 The fact that higher-order differential operators depending on e may take a simpler form after e-expansion has also been 
observed in the evaluation of multi-loop master integrals; see e.g. Refs. [46, 48-50, 53. |62H 
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and lfc = 1. As a first step toward the construction of a solution, the operator should be rewritten in the 
form 



(p) _ 



where ^ are rational numbers. Eqs. (|55|) and (|36l) give rise to a system of equations 



b*p$~hp) + p r(p) = pr L (B-i), k = i, • • • , P . 



R 1 pc 1 1 (P)+pr l (^ = p p k (A), k=i, 



The differential equation 



generates the function 



-(l-z) —+Ri-^R 2 } //( :! =0. 



h(z) = Cz- R *(z-l) R2 - Ri , 
for which only three rational parametrizations are known: 

(i)R 1 = R 2 ; {ii)Ri=0; (m)R 2 = 



(36) 



(37) 



(38) 



(39) 



(40) 



When all A; = and Bj = 1, we have integer parameters and the differential operators are factorizable. 
When only one upper parameter A\ ^ and one lower parameter B\ ^ 1, the f3j again all vanish, and a 
parametrization z — > £(z) should exist such that [(1 — z)h(z)]~ 1 dz/ d£., [zh(z)]~ 1 dz/d£ t , z~ 1 dz/d£ t , are rational 
functions of £, where /i(z) = Cz 1 " 51 ^ - (See Ref. [Hjj.) 



Example III: Let us consider a Gauss hypergeometric function of the form uj(z) = 2-Fi( — + ai£, — + a2£; 
1 — - + ce; z) , where p\,p 2 , q, T are integers. It is a solution of the differential equation 



d pi 

z— H h aie 

az q 



z- — | h a 2 e I ui(z) 



dz q 



r 



dz \ dz 



ce I w(z) 



(41) 



for coefficient functions U!k(z) defined via the expansion oj(z) — 1 + X)^Li ^fcOO^ • Eq. (1411) produces an infinite 
system of linear differential equations 



Pi +P2 



- 0+~ 



Pi 



q 



P2 



q 



i 



9 



Wfc(z) 



= (ai + a 2 -- ) ( ) Wfe-i(z) + -fiuk-i 



P2 



Pi 



Ldk-i(z)+aia 2 UJk-2(z), (42) 



where (3 is arbitrary and k runs from to oo. The second order differential equation (|42|) can be split into two 
first-order differential equations by introducing new functions 



Pk{z) = ( z -^ + P ) > 



(43) 



so that 



t^ \ d , fa Pi +P2 



PfeW 



P- 



Pi 



P2 



1 



CJk(z) 



= (ai+a 2 -^j pk-i(z) + ^fiuik-i - ai\B - u> k -i{z) - a,2 \J3 - ^-J ijj k -i{z)+aia 2 uJk-2{z). (44) 

We can use the fact that this system takes triangular form when the last term in the l.h.s. of the second equation 
is zero to obtain the following solutions: 



P1P2 =0-^/3 = 0, 
Pl = — > = , 

L - El 
q q 



p 2 =p 



(45) 
(46) 

(47) 
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We would like to analyze the case of the parameters in Eq. (|47j) in detail. Then 

p k {z) = (z^--- ) u) k (z) , (48) 



dz q 
d p 



dz q 



(c\ c p\ J)\ — p 
a 1 +a 2 ~-j pk-i{z)-\ w fc _i(z)-ai Wk-\{z)+a 1 a 2 U3 k -2{z) ■ (49) 
z / z c[ q 

The redefinition (wfe, p^) — > (^z^ ir k , (1 — z)~ « 07^ leads to a new set of equations 

(7 fe (^) = h(z)z-^-n k (z) , (50) 



(l-z)-^-cr fc (2;) = fai+a 2 -- N ) (T fe _i(^) + -— h(z)n k -i(z)-ai— -h(z)i: k ^ 1 {z)+a 1 a 2 h(z)'K k ^ 2 (z),{bl) 

dz \ z) z q q 



where h(z) = (1 - z)P^z r / c '. The result can be expressed in terms of Goncharov's polylogarithms if there is a 

parametrization z — > £(z) such that [zh(z)]~ 1 dz/d^, [(1 — z)h(z)]~ 1 dz/d^, z~ 1 dz/d£, and (1 — z)~ 1 dz/d£ > are 

■ — i ... / \Va 

rational functions. (See Ref. |24|.) Such a parametrization exists when p = — r and z — > £ = I — ) . In this 

way, we find the following set of conditions: 

p = _r_ = Pl = P2 j 
9 9 9 

Under these conditions, the coefficients of the £-expansion of the function 2 Fi(^r + o,\e, ^ + a 2 e; 1 — ~ + ce; z) , 
are expressible in term of Goncharov's polylogarithms. Eq. (|5"!2)) corresponds to Lemma IV of Ref. [23J, but the 
present derivation does not rely on any symmetries of Gauss hypergeometric functions. 

Remark 1: For a Gauss hypergeometric function, the e = term should be a rational (not just algebraic) 
function. 

Remark 2: For Eq. (l4"B"l) . an additional relation between p 2 and r arises from Eq. ([58)1 . and we arrive again at 
one of the cases of Eq. (l4*Ul) . 

Let us consider the e expansion of a hypergeometric function with the following set of parameters: uj(z) — 
pF p -i (as, Ai + cie, A 2 +c 2 s; be, Bi+f±e, B 2 +f 2 s; z), where A±,A 2 , B\,B 2 , a, b, c, and / are rational numbers. 
Eqs. ([37| take the form 

R 2 + (3 = (fli - 1) + (B 2 - 1) , R t + = A x + A 2 , 

R2P = (Si - l)(B 2 - 1) , RiP = A X A 2 . (53) 

Eq. (|4"0")l provides additional conditions on the relations between R\ and R2, and as a consequence, there are 
three different solutions: 

(54) 
(55) 
(56) 

The solutions of Eqs. ([5~4")l and ([5"6")l are the roots of a quadratic equation — (A\ + A 2 )x + A1A2 = 0, and the 
solution of Eq. (I55[) satisfies the same quadratic equation with the replacement Ai — ?> Bi — 1. One solution of 
Eq. ([51]) is ; = 1 + Aj, and one solution of Eqs. ([55} and ((56j) is B x = 1 + A x . 

(p) 

There is another parametrization for the same hypergeometric function. Let us rewrite the operator Dq in 



Ri — R2 ■ 




- 1 + 


B 2 ~ 


l = ^i 


+ A 2 , (B 1 -l)(B 2 -l)=AiA a ; 


Rt = : 


A 2 


= , 


P = 


Ai , 






Ri 


= (Bi 


-1) 


+ (£ 2 - 


1)-Ax, fl 2 ^i = (Bi - 1){B 2 - 1 


R2 =0 : 


B, 


= 1 , 


= 








Pi 


= A X 


+ A 2 




1), i? 1 (B 1 -l)=A 1 A 2 ; 
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Eq. (SD as 



fe=i 



•A 1 A 2 -^(Bi-l)(B2-l)>e p - 2 



-(l-z)^-(i8-yli-A 2 )+i {p-B x -B 2 +2) 
1 



(0+/3) 



,p-2 



(57) 



+ \{fi- A x ) (13 -A 2 )--(f3- G8-.B a +l) }> # p - 2 



The first condition is that there should exist a common factor for the last line in Eq. (1571) . for example, 
B\ = A\ + 1, and consequently, /3 = Ai. 



Example IV: Let us consider the e expansion of a hypergeometric function 3F2 with the following set of param- 
eters: lj(z) — 3F2 (^ + aie, C12E, a 3 e; l + ^+&i£, 1 — 1+&2£; z^j with 01,02,03 7^ 0. Starting from the differential 
equation for this hypergeometric function, 



z ft 






>+-+ai£ ) 




9 / 



p 



w(z) = . 



(58) 



the system of differential equations for the coefficient functions Wk{z) in its e expansion ui(z) = l+^Z^i w k(z)e k 
is found to be 



(9 + - ) 0ui m (z) = 



h + b 2 



r \ r 
+ - J 6»w m _i(z) + a 2 a 3 ~cj m _ 2 (z) 



+aia 2 a 3 w m _ 3 (z) + <5i0w m _i(,z) + 5 2 0uj m -2(z) + -<5 3 6»a; m _i(z) + -&t0u; TO _ 2 (z) 

z z 



where 8 — zd/dz and Sj are constants. After a redefinition 



cj fc (z),6>cj fc (z), ( 6 + - ) 9uj k (z), 



uj k (z),z r/q a k (z), 



z - 1 



p/g 



■Az) , 



(59) 



(60) 



we find that the function h(z) defined by (f3"9")l is h(z) = z^ p+r ^ q (z~ \)~ p / q . By (|4T))) . a rational parametrization 
is possible when p = — r. 

The methods described here can be extended to any multi-loop Horn-type hypergeometric function. Starting 
from the Pfaff form of the differential equation, 



*W (z, e) = Y, PiJ,k (z, e)^ (z, e)dz k 



(61) 



where Pij^ k (z,s) are rational functions, the system can be transformed to triangular form and integrated. 

Let us consider the e-expansion of the Appell hypergeometric function F3, which was analyzed in the context 
of photon box diagrams |63j : 

f Pi P2 T\ f2 P 
F 3 hOie, ha 2 e, \-o\E, ho 2 e, 1 \- ce;x,y 



q q 

00 00 

EE 

m=0 n— 



x m y n 



,1 «! 



(62) 



m+n 



Applying our methods, we find that the coefficients of e-expansion may be expressed in terms of Goncharov's 
poly logarithms when pjrj = for j = 1,2, and a rational parametrization should exist for the functions 



(x - l) s i+P 

H{x,y) = (-l)( s i+^)Az 



1/9 



X S2+PySl+P 



h»(x) = (-iy^ q 
1/9 



(y - i) S2+p 



1/'/ 



(xy-X- y) s i+ s 2+P 



(63) 
(64) 
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where Sj = pj + Tj and j = 1, 2. As result of our analysis, we claim that in two cases (only), 

F 3 [ h + — + ai£,I 2 +a 2 e,I 3 + bi£,h + b 2 e,I 5 + — +ce;x,y ) , (65) 

\ q q J 

F s (h + ^- + h + a 2 £, h + he, h + b 2 £, h + ce; x, y\ , (66) 

with integer values Ij,pi,q, the e expansion can be expressed in terms of Goncharov polylogarithms. We 
note that for F3, the e = term should be a rational function. 
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